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I present two recent pieces of work on semileptonic B ii, p decays where it is shown how lattice QCD 
calculations can be used to test heavy-quark symmetry and determine phenomenologically relevant quantities 
despite the limits on these calculations' kinematical reach. The study of semileptonic B —> p decays was performed 
with the UKQCD Collaboration. 



1. Motivations 



spin effects can be measured. 



The CLEO Collaboration has very recently pre- 
sented measurements of the branching ratios for 
B —^ piv and B —> Trii? decays {i=e, fj,) |^ and 
has announced its intention to measure the corre- 
sponding differential decay rates. These various 
measurements represent an excellent opportunity 
to determine the poorly known CKM matrix el- 
ement \Vub\- Such determinations require under- 
standing the non-perturbative, strong-interaction 
corrections to the elementary b — u — W coupling 
contained in the matrix elements of the weak cur- 
rents V^=uj'^b and A'^=u'-f^j^b between B and 
TT or p meson states. It is to calculate these matrix 
elements that we resort to the lattice. 

heavy light quark decays, such as the ones 
that concern us here, are also interesting because 
they enable one to test heavy-quark symmetry 
(HQS). For these decays, HQS is weaker than 
for heavy — > heavy quark decays: it only ap- 
plies in a limited region around the zero recoil 
point q'^=(lmaxi where q is the four-momentum 
transferred to the leptons, and imposes no nor- 
malization condition on the relevant form factors 
at Qmax- Nevertheless, because both the mass and 
the spin of the heavy quark can be varied in lat- 
tice calculations, the deviations from the heavy- 
quark limit due to finite heavy-quark mass and 



2. Limitations 
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Current day lattice calculations. 



with lattice 

spacings on the order of 3 GeV^^, do not per- 
mit one to simulate b light quark decays over 
their full kinematical range. The problem is that 
the energies and momenta of the particles in- 
volved, whose orders of magnitude are set by the b 
quark mass (nib ~ 5 GeV), are large on the scale 
of the cutoff in much of phase space. To limit 
these energies in relativsitc lattice quark calcu- 
lations, one performs the simulation with heavy- 
quark mass values mg around that of the charm 
(rric ~ 1.5 GeV), where discretization errors re- 
main under control. Then one extrapolates the 
results up to mi, by fitting heavy-quark scaling 
relations (HQSR) with power corrections to the 
lattice results (see Sec. 3.1). Another approach is 



to work with discretized versions of effective the- 
ories such as Non-Relativistic QCD (NRQCD) or 
Heavy-Quark Effective Theory (HQET) in which 
the mass of the heavy quark is factored out of 
the dynamics. All approaches, however, are con- 
strained to relatively small momentum transfers 
because of the limited applicability of HQS and 
because of momentum-dependent discretization 
errors. So one can only reconstruct the de- 
pendence of the relevant form factors in a limited 
region around q'^ax ^^'^ '^^^ is left with the prob- 
lem of extrapolating these results to smaller q^ . 

heavy — > light quark decays are difficult 
in any theoretical approach. Indeed, they re- 
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quire understanding the underlying QCD dynam- 
ics over a large range of momentum transfers 
from g^„^=26.4 GeV^(20.3 GeV^) for semilep- 
tonic B — > TT (i? ^ p) decays, where the final 
state hadron is at rest in the frame of the B me- 
son, to where it recoils very strongly. 

3. 5° p^i^D and a Model-Independent 
Determination of \Vub\ o 

One solution to the problem of the limited kine- 
matical range of the lattice results is to ignore 
the problem or rather rely on the ingenuity of 
experimental groups to provide measurements of 
partial rates in the region where lattice results 
are available. Combined with lattice results for 
B — !■ piD decays, such experimental measure- 
ments will enable a model-independent determi- 
nation of \Vub\ [§]■ Rates should be sufficient 
since our lattice results span a range of from 
~ 14.4 GeV^ to ql„^ over which the partially in- 
tegrated is 4.6 t3|Ki(,pps^^. This represents ap- 
proximatively 1/3 of the total rate obtained from 
light-cone sumrules (LCSR) in Ref . |Q , whose re- 
sults at large agree well with ours. 

3.1. Form Factors and Heavy-Quark Ex- 
trapolation 

To describe B^ p^£~D decays, we must 
evaluate the matrix element {p~^ {p' , r])\V — 
A'^\B°{p)), traditionally decomposed in terms of 
four form factors Ai, A2, A and V which are 
functions of q"^, where q=p — p' . We calculate 
this matrix element for four values of the heavy- 
quark mass around that of the charm. Then, to 
obtain A\ at the scale of the B meson we fit, to 
the lattice results, the HQSR 

Ai(u;,M)a,(M)2/'3oVM =cH(l + %i 

+0(A2/Af2)), (1) 

where M is the mass of the decaying meson, 
/5o=ll — Irifl?) and A is an energy characteris- 
tic of the light degrees of freedom. This scaling 
relation holds for a;=(M^ + — q^) /2Mm close 

^The results presented here were obtained on a 24'^ X 48 
lattice at /3 = 6.2 from 60 quenched configurations, using 
an O (a)-improved SW action for the quarks. (Please see 
for details.) 
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Figure 1. Ai vs. q^ from UKQCD (crosses), 
APEg (right-hand dot), ELC[| (left-hand dot) 
andLCSR[§ (curve). (Adapted from Ref. §.) 



to 1 and the fit parameters c and d are indepen- 
dent of M. Here m is the mass of the final state 
meson. Once c and d are fixed, it is trivial to 
obtain Ai{M=mB) at the corresponding value of 
CO. Furthermore, d determines the size of correc- 
tions to the heavy-quark limit. Repeating this 
procedure for many values of lu, one obtains the 
q^ dependence of the desired form factor around 
Imax- The resulting Ai{q^,mB) is plotted to- 
gether with the LCSR result of Ref. P] and the 
lattice results of APE and ELC Agree- 
ment amongst these four calculations is excellent 
as it for A2 and V [||J^ , which are obtained in an 
entirely analogous way. 

3.2. Rates 

Having determined Ai , A2 and we can com- 
pute {1 / \ Vub\'^)dT / dq'^ . Our results are plotted in 
Fig. H (squares). In the region of q^ accessed, we 
can legitimately expand, around f/^^j., the helic- 
ity amplitudes that appear in the rate. Thus we 
fit, to the lattice points, the parametrization 

^Q'" ~ 91 2^1/2 

X a'{l + b{q'-ql,J) (2) 

where \{q^) = {]VP + w? - q^) - AAPm^. We find 
a=4.6 tt ± 0.6 GeV and 6=(-8 ^1)1^-"^ GeV"^ 
where the second error on a is systematic, all 
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Figure 2. The data points are our lattice results 
and the solid curve, the fit to Eq. (0). 



other errors being statistical. With a and b de- 
termined, the only unknown in Eq. (|^) is \Vub\- 
Therefore, a fit of the parametrization of Eq. (|^) 
to an experimental measurement of the differen- 
tial decay rate around q^ax determines \Vub\- In 
this determination, a plays the role of in 
the extraction of |Vc&| from semileptonic B ^ D 
or D* decays [§ and b the role of The 
difference, here, is that a is not determined by 
HQS up to small radiative and power correc- 
tions. It is a genuinely non-perturbative quan- 
tity. Another way of determining \Vub\ from the 
lattice results is to compare partially integrated 
rates from > 14 GeV^ to g^^^ given by Eq. (||) 
to the corresponding experimental measurements. 
Both these methods yield \Vub\ with approxima- 
tively 10% statistical and 12% theoretical uncer- 
tainties. 

3.3. A Test of HQS 

In Fig. ^ we compare semileptonic B ^ p 
form factors with those governing the short dis- 
tance contribution to radiative B K*j de- 
cays for which the relevant hadronic matrix ele- 
ment is {K*{p\r])\saf"'q''bR\B{p)), with q=p-p'. 
This matrix element is traditionally decomposed 
in terms of three form factors, Ti, T2 and T3. 
The comparison is made for three initial meson 
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Figure 3. Ratios V/2Ti and Ai/2iT2 for 5 values 
of to and three initial meson masses. The solid 
lines are the HQS predictions. 



masses: AI—m£), M—ms and M — s- 00. For 
identical final-state vector mesons (in Fig. I aU 
light-quarks involved have the same mass, slightly 
larger than that of the strange), HQS predicts 



V{q^) = 2Ti(q2), A, = 2tT2{q^) , 



(3) 



for q^ around g^^j^. or, equivalently, to close to 1. 
While V/2Ti displays large 1/M corrections at 
the D and even B meson scale, Ai/2iT2 exhibits 
no such corrections even at the D scale. Both 
ratios, however, converge to 1 in the heavy-quark 
limit which gives us confidence that we control 
the heavy-quark-mass dependence of the various 
form factors. Furthermore, these ratios can help 
constrain the possible q^ dependences of the var- 
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ious form factors around q^^ax M=mB- 

4. 5° TT^t^v and Dispersive Constraints 

A second solution to the problem of the lim- 
ited kinematical reach of lattice simulations of 
heavy light quark decays is to combine lat- 
tice results for the relevant form factors around 
^max with dispersive bound techniques to ob- 
tain improved, model-independent bounds for the 
form factors for all For the case of 

& TT^£~D decays, whose hadronic matrix ele- 
ment, {tt^ {p')\V^^\B° (p)) , is traditionally decom- 
posed in terms of two form factors f^{q^) and 
f^{q^): one can use the kinematical constraint, 
/+(0)=/°(0), to further constrain the bounds. 



4.1. Dispersive Bounds 

The subject of dispersive bounds in semilep- 
tonic decays has a long history going back to 
S. Okubo et al. who applied them to scmilcp- 
tonic K t: decays Q. C. Bourrely et al. first 
combined these techniques with QCD and ap- 
plied them to semileptonic D ^ K decays [10 . 
Very recently, C.G. Boyd et al. applied them to 
B -Kiv decays 0. 

The starting point for B tt£D decays is the 
polarization function 



Il^'''{q)^^ J d^x e^9-^(o|r(t/^(a;)l/''t(o)) |o) 

Hl^q" - g'-^'q^) Tlriq^) + q'^q'' n^iq^) , (4) 

where IVt(l) corresponds to the propagation of a 
J^=l^ (0+) particle. The corresponding spectral 
functions, 1u\T\t,Li are sums of positive contribu- 
tions coming from intermediate B* ( J^=l~), Bit 
[J^^Q^ and 1^), ... states and are thus upper 
bounds on the Bt: contributions. Combining, for 
instance, the bound from Im 11^ with the disper- 
sion relation {Q^— — q^) 



aQ2 



tluvULit) 

{t + Q2)2 



dt 



one finds 



1 /■°° 

XLiQ')>- dtk{t,Q')\f{t)\\ 



(5) 



(6) 



where t± = {mB^'rnT^)'^ and k{t,Q^) is a kinemat- 
ical factor. Now, since Xl{Q^) can be calculated 
analytically in QCD for far enough below the 
resonance region (i.e. —Q^ <C m^), Eq. (||) gives 
an upper bound on the weighted integral of the 
magnitude squared of the form factor /° along the 
Btt cut. To translate this bound into a bound on 
/" in the region of physical B iriD decays is a 
problem in complex analysis (please see Ref. Q 
for details). A similar constraint can be obtained 
from IIt- for /+. There, however, one has to con- 
front the additional difficulty that /+ is not an- 
alytic below the Bn threshold because of the B* 
pole. 

The beauty of the methods of Ref. ||T^ is that 
they enable one to incorporate information about 
the form factors, such as their values at various 
kinematical points, to constrain the bounds. For 
the case at hand, however, these methods must be 
generalized in two non-trivial ways. In construct- 
ing these generalizations, one must keep in mind 
that the bounds: 1) form inseparable pairs; 2) do 
not indicate the probability that the form factor 
will take on any particular value within them. 

4.2. Imposing the Kinematical Constraint 

The first problem is that Eq. (^ and the equiv- 
alent constraint for /+ yield independent bounds 
on the form factors which do not satisfy the kine- 
matical constraint /~''(0)=/°(0). The bounds on 
require /^(O) to lie within an interval of val- 
ues /+ and those on within an interval Iq. To- 
gether with these bounds, however, the kinemat- 
ical constraint requires /'*'(0)=/''(0) to lie some- 
where within /_!_ n /q. Thus, we seek bounds on 
the form factors which are consistent with this 
new constraint. 

A natural definition is to require these new 
bounds to be the envelope of the set of pairs of 
bounds obtained by allowing f~^{0) and /"(O) to 
take all possible values within the interval I+DIq. 
In Ref. |8| , I show how this envelope can be con- 
structed efficiently and that the additional con- 
straint can only improve the bounds on the form 
factors for all q^. Also, as a by product, one ob- 
tains a formalism which enables one to constrain 
bounds on a form factor with the knowledge that 
it must lie within an interval of values at one or 
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more values of . 

4.3. Taking Errors into Account 

As they stand, the methods of Ref. can only 
accommodate exact values of the form factors at 
given kinematical points and contain no provi- 
sions for taking errors on these values into ac- 
count. Of course, the results given by the lattice 
do carry error bars. More precisely, the lattice 
provides a probability distribution for the value 
of the form factors at various kinematical points. 
What must be done, then, is to translate this dis- 
tribution into some sort of probability statement 
on the bounds. The conservative solution is to 
consider the probability that complete pairs of 
bounds lie within a given finite interval at each 
value of . Then, using this new probability, one 
can define upper and lower p% bounds at each q^ 
as the upper and lower boundaries of the interval 
that contains the central p% of this probability.H 
These bounds indicate that there is at least a p% 
probability that the form factors lie within them 
at each q^ . 

4.4. Lattice-Constrained Bounds 

To constrain the bounds on /+ and /°, I use the 
lattice results of the UKQCD Collaboration [|l2), 
to which I add a large range of systematic errors 
to ensure that the bounds obtained are conser- 
vative. Because of these systematic errors, the 
probability distribution of the lattice results is 
not known. I make the simplifying and rather 
conservative assumption that the results are un- 
correlated and gaussian distributed. I construct 
the required probability by generating 4000 pairs 
of bounds from a Monte-Carlo on the distribution 
of the lattice results. My results for the bounds 
on the form factors are shown in Fig. ^. I have 
plotted the two form factors back-to-back to show 
the effect of the kinematical constraint. With- 
out it, the bounds on /+ would be looser, espe- 
cially around g^=0, where phase space is large. 
Since /"*" determines the rate, the kinematical 
constraint and the bounds on /° are important. 

Also shown in Fig. ^ is the LCSR result of 

■^The density of pairs of bounds increases toward the center 
of the "distribution" as long as the distribution of the 
lattice results does. 




Figure 4. and f^{q^) versus q^ . The 

data points are the lattice results of UKQCD jl2j 
with added systematic errors. The pairs of fine 
curves are, from the outermost to the innermost, 
the 95%, 70% and 30% bounds. The shaded curve 
is the LCSR result of Ref. 



Ref. |13[ which has two components: for q^ below 
15 GeV , the q^ dependence of is determined 
directly from the sumrule; for larger q^ , pole dom- 
inance is assumed with a residue determined from 
the same correlator. Agreement with the bounds 
is excellent. In Ref. the bounds are compared 
with the predictions of more authors as well as 
with direct fits of various parametrizations to the 
lattice results. Though certain predictions are 
strongly disfavored, the lattice results and bounds 
will have to improve before a firm conclusion can 
be drawn as to the precise q^ dependence of the 
form factors. 

The bounds on /+ also enable one to con- 
strain the B*B7r coupling gs-BTT which deter- 
mines the residue of the B* pole contribution to 
/"*". The constraints obtained are poor because 
/+ is weakly bounded at large q^ . Fitting the 
lattice results for /" and /+ to a parametriza- 
tion which assumes B* pole dominance for 
and which is consistent with HQS and the kine- 
matical constraint gives the more precise result 
3B*+B°Tr+ = 28 ± 4.B However, because this re- 
sult is model-dependent, it should be taken with 

•^The result of this fit is entirely compatible with our 
bounds on /+ and 
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Table 1 

Bounds on rate in units of iKibpps^-'^ 
/+(0). 



and 
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95% 


2.8 - 


-> 24 


-0.18 


-> 0.85 


90% 


3.6 - 
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0.00- 


0.68 


70% 


4.4- 


13 


0.10 - 


0.57 


50% 


4.8 - 


10 


0.18 - 


0.49 


30% 



4.5. Bounds on the rate and \Vub\ 

As was done for the form factors, one can de- 
fine the probabihty of finding a complete pair of 
bounds on the rate in a given interval and from 
that probability determine confidence level (CL) 
intervals for the rate. The resulting bounds are 
summarized in Table |l|. They were obtained by 
appropriately integrating the 4000 bounds gener- 
ated for f^iq^), taking the skewness of the re- 
sulting "distribution" of bounds on the rate into 
account. The CL bounds obtained can be used, 
in conjunction with the branching ratio measure- 
ment of CLEO Ipl, to determine \Vub\- One finds 



|Kf,|10V^B°/l-56ps- (34-49)±8±6 , (7) 

where the range given in parentheses is that ob- 
tained from the 30% CL bounds on the rate and 
represents the most probable range of values for 
\Vub\- The first set of errors is obtained from 
the 70% CL bounds and the second is obtained 
by combining all experimental uncertainties in 
quadrature and applying them to the average 
value of \Vub\ given by the 30% CL results. This 
determination of \Vub\ has a theoretical error of 
approximately 37%. Though non-negligible, this 
error is quite reasonable given that the bounds on 
the rate are completely model-independent and 
are obtained from lattice data which lie in a lim- 
ited kinematical domain and include a conserva- 
tive range of systematic errors. 

5. Conclusion and Outlook 

Because HQS applies to heavy — > light quark 
decays in a rather limited way, it is not possible 



to determine the full dependence of the rele- 
vant form factors from the lattice alone. The fiip 
side of the coin is that the model-independent in- 
formation provided by lattice calculations about 
these decays, though limited, is still very impor- 
tant, because the relevant matrix elements are 
not anchored at zero recoil by HQS, up to small 
radiative and power corrections, as they are in 
heavy — *■ heavy quark decays. 

I have presented two approaches by which the 
information provided by the lattice on exclu- 
sive semileptonic b ^ u decays can be used 
to extract \Vub\- Both approaches will bene- 
fit from forthcoming, improved lattice results. 
The lattice-constrained bounds would also ben- 
efit enormously from an increase in the range of 
accessible q^. 

Finally, the techniques developed in Ref. Q to 
construct lattice-improved bounds for B — > niD 
decays are in principle applicable to limited re- 
sults obtained by non-lattice means and to other 
processes such as i? ^ p£D and B K*^ decays. 

REFERENCES 



[1] 
[2] 

[3] 



[4] 
[5] 
[6] 



[7] 
M. 



CLEO Collaboration, CLEO CONF 96-16 
UKQCD Collaboration, J. M. Flynn et i 
Nucl. Phys. B461 (1996) 327. 
P. Ball, talk given at the XXXIst Rencontres 
de Moriond, "Electroweak Interactions ", hep 
ph/9605233| . 



APE, Phys. Lett. B345 (1995) 513. 
ELC, Nucl. Phys. B416 (1994) 675. 
J.M. Flynn, "Heavy Light Weak Matrix Ele- 
ments" in Proc. 14th Int. Symp. on Latt. Field 
Theory, St Louis, USA, 4-8 June 1996. 

M. Neubert, Phys. Lett. B338 (1994) 84^ 

L. Lellouch CPT-95/P.3236 (revised), |ep| 
ph/9509358| , to appear in Nucl. Phys. B. 
[9] S. Okubo and I-Fu Shih, Phys. Rev. D4 
(1971) 2020. 

[10] C. Bourrely, B. Machet and E. de Rafael, 

Nucl. Phys. B189 (1981) 157. 
[11]C.G. Boyd, B. Grinstein and R.F. Lebed, 

Phys. Rev. Lett. 74 (1995) 4603. 
[12] UKQCD, Nucl. Phys. B447 (1995) 425. 
[13] V.M. Belyaev, V.M. Braun, A. Khodjamirian 

and R. Riickl, Phys. Rev. D51 (1995) 6177. 



